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A  REGION  OF  SATISFACTORY  PERIOD-DAMPING  RELATIONSHIP 
OF  THE  OSCILLATORY  MODE  OF  MOTION 
By  Leonard  Sternfield  and  Ordway  B.  Gates,  Jr. 

SUMMARY 


A  method  has  Been  derived  hy  which  a  boundary  can  be  obtained  that 
defines  a  region  in  which  there  exists  a  satisfactory  relationship 
between  the  period  and  damping  of  the  lateral  oscillatory  mode  of  motion, 
according  to  any  given  criterion  for  this  relationship.  In  addition, 
a  method  is  discussed  by  which  curves  representing  a  constant  rate  of 
spiral  divergence  may  be  constructed. 

The  methods  aB  presented  are  applicable  to  both  lateral-stability 
and  longitudinal -stability  analyses. 


INTRODUCTION 


In  lateral-stability  analyses  it  is  the  usual  practice  to  calculate 
a  neutral-oscillatory-stability  boundary  that  is  plotted  as  a  function 
of  the  directional-stability  derivative  Cnp  and  the  effective-dihedral 

derivative  C7  .  This  boundary  has  a  special  significance,  since  for  a 
*0 

particular  airplane  it  indicates  the  combinations  of  C  and  C, 

n0  Z0 

necessary  for  oscillatory  stability.  Although  this  boundary  is  definitely 
an  aid  in  a  stability  analysis,  its  value  is  somewhat  restricted  in  that 
it  affords  no  information  about  the  variation  of  the  period-damping 
relationship  throughout  the  stable  region.  That  is,  that  region  of 
the  C  ,C7  plane  in  which  a  given  criterion  for  satisfactory  relationship 
3  6  3 

of  the  period  and  damping  of  the  oscillatory  mode  will  be  satisfied 
cannot  be  determined  from  this  boundary.  This  information  can  be  obtained 
by  using  the  methods  described  in  reference  1  to  calculate  curves  of 
constant  period  and  constant  damping,  but  because  of  the  complexity  of 
the  calculations  involved,  a  more  direct  approach  to  the  problem  is 
desirable.  In  the  present  paper  a  method  is  derived  by  which  it  is  possible 
to  obtain  directly  a  boundary,  plotted  as  a  function  of  C  and  C-, 

Hd  la  9 
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that  will  define  a  region  in  which  there  exists  a  satisfactory  relation¬ 
ship  between  the  period  and  damping  of  the  oscillatory  mode  of  motion, 
according  to  a  prescribed  criterion  for  this  relationship. 

For  the  combinations  of  CL  and  C,  necessary  to  satisfy  various 

criterions  for  the  period  and  damping  of  the  oscillatory  mode,  undesirable 
spiral  instability  may  be  present.  In  view  of  this  fact  a  method  is  also 
presented  by  which  curves  representing  a  constant  rate  of  spiral  divergence 
may  be  constructed. 

SYMBOI£>  AND  COEFFICIENTS 
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angle  of  roll,  radians 
angle  of  yaw,  radians 
angle  of  sideslip,  radians  (v/v) 

sideslip  velocity  along  the  Y-axis,  feet  per  second 
airspeed,  feet  per  second 
mass  density  of  air,  slugs  per  cubic  foot 
dynamic  pressure,  pounds  per  square  foot 
wing  span,  feet 
wing  area,  square  feet 
weight  of  airplane,  pounds 
mass  of  airplane,  slugs  (w/g) 
acceleration  due  to  gravity,  feet  per  second  per  second 
relative -density  factor  (m/pSh) 

inclination  of  principal  longitudinal  axis  of  airplane  with  respect 
to  flight  path,  positive  when  principal  axis  is  above  flight 
path  at  the  nose,  degrees 

angle  of  flight  path  to  horizontal  axis,  positive  in  a  climb,  degrees 
radius  of  gyration  in  roll  about  principal  longitudinal  axis,  feet 
radius  of  gyration  in  yaw  about  principal  vertical  axis,  feet 
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K, 


K, 


% 


K- 


XZ 


CL 

cz 

Cn 


'n, 


'Yp 


'XI, 


'np 


nondimensional  radius  of  gyration  in  roll  about  principal 
longitudinal  axis  ^kx^/b^ 

nondimenaional  radius  of  gyration  in  yaw  about  principal 
vertical  axis 

nondimensional  radius  of  gyration  in  roll  about  longitudinal 
stability  axis  ^ 2cos2ti  +  2sin2T]  ^ 

nondimensional  radius  of  gyration  in  yaw  about  vertical 
stability  axis  (K  ^cos^r)  +  Kj  ^sin^T]  ^ 

nondimensional  product-of -inertia  parameter 

((Kz02  ‘  Kx02)  sin  *  003  T>) 

trim  lift  coefficient  f—  003  2L\ 

V  <1S  ) 

rolling -moment  coefficient  ( ^ingmomentN 

\  <lSb  ) 

yawing -moment  coefficient  ^  Yawlng^ome  nt  ^ 

lateral-force  coefficient  (  Jia^era-1..  Yorce \ 

\  q.S  J 

effective -dihedral  derivative,  rate  of  change  of  rolling -moment 
coefficient  with  angle  of  sideslip,  per  radian  (dC^dp) 

directional-stability  derivative,  rate  of  change  of  yawing -moment 
coefficient  with  angle  of  sideslip,  per  radian  (dC^dp) 

lateral-force  derivative,  rate  of  change  of  lateral-force 
coefficient  with  angle  of  sideslip,  per  radian  (dCy^dp) 

damping -in -yaw  derivative,  rate  of  change  of  yawing-moment 
coefficient  with  yawing -angular -velocity  factor,  per 


radian 


rV 


■orb 


rate  of  change  of  yawing-moment  coefficient  with  rolling -angular - 


velocity  factor,  per  radian 


damping-in-roll  derivative,  rate  of  change  of  rolling -moment 
coefficient  with  rolling-angular-velocity  factor,  per 


radian  I  dC 


>/*§) 


k 


NACA  TN  No.  1859 


rate  of  change  of  rolling  "moment  coefficient  with  yawing -angular 
velocity  factor,  per  radian 


rate  of  change  of  lateral-force  coefficient  with  rolling -angular 
velocity  factor,  per  radian 


rate  of  change  of  lateral-force  coefficient  with  yawing -angular - 
velocity  factor,  per  radian 

rate  of  change  of  rolling -moment  coefficient  with  angle  of  roll, 
per  radian 

time,  seconds 

nondimens ional  time  parameter  based  on  span  (Vt/b) 


differential  operator 
Routh’s  discriminant 


complex  root  of  stability  equation 


AX**  +  B\3  +  CA,2  +  DX  +  E  =  0  l  X±  =  a  +  ico  =  Re 


R  =  va^  + 


Q  =  tan' 


factors  used  in  calculation  of  boundary  of  satisfactory  period- 
damping  relationship  fe-n  =  CQS  for  0  <  n  ^  k\ 


damping  relationship  Kn  =  CQ^nn^  $n 


sin  n 9 
Rn-1 


X2  spiral  stability  root  of  stability  equation 

AxA  +  BX3  +  CX2  +  DX.  +  E  =  0 
P  period  of  oscillation,  seconds 

Tjl j2  time  for  amplitude  of  oscillatory  or  spiral  mode  to  decrease  to 
one -half  its  original  value 

T  time  for  amplitude  of  oscillatory  or  spiral  mode  to  increase  to 

L  double  its  original  value 


A,B,C,D,E  coefficients  of  lateral -stability  equation 
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EQUATIONS  OF  MOTION 


The  nondimens ional  linearized  equations  of  motion,  referred  to  the 
stability  axes,  used  to  calculate  the  spiral- stability  and  oscillatory- 
stability  boundaries  for  any  flight  condition,  are: 

Rolling 

+  KxzDb2t^=  +  ^CzrUb^  +  ^CipDb0 

Yawing 

^b^zV*  +  %zVV)  =  CnpP  + 

Sideslipping 

2Hb^P  +  V)  =  °YpP  +  +  CL$  +  1%^ *  +  (CL  tan  7^ 

When  ls  substituted  for  $ ,  'foe^8^  for  i,  and  P0e^3^  for 

in  the  equations  "written  in  determinant  form,  X  must  be  a  root  of  the 
stability  equation 


AX1*  +  BX3  +  CX2  +  DX  +  E  =  0 


(1) 


where 


A  =  8nb3(Kx2Kz2  -  KXZ2) 

B  =  -2pt2  (2Kx2Kz2CYp  +  Kx2Cnr  +  -  2Kxz2CY(3  - 

c  =  »ib(K:X2cnrcYp  +  ^bKX2°np  +  KZ2c2pCYp  +  |CnrC  Zp  "  KXZCZrCYp 

-  V"b%zClp  -  CnpKxzCYp  -  ^npClr  +  KXZCnpCYp  "  KZ2cYpClp 

-  KX2CyrCnp  +  KxZcYrC2p) 
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D  =  -tCnrcZpcYp  -  ^CZpCnp  +  jj?npCZrcYp  +  ^bCnpCZp  +  ^\>CLKJZ% 


-  2nbCLKz2cZp  -  2^Kx2CnpcL  tan  7  +  2nbKxzCzpCL  tan  7 
+  -  JCnpCZpCYr  “  t°ZrC npCYp  +  £%%% 

E  =  2CL(Cnr%  "  GzrCnp)  +  2°L  7(CzpCn3  "  %%) 

The  damping  and  period  of*  the  lateral  oscillation  in  seconds  are  given 
"by  the  equations 


Tl/2 


-0-693  b 
a  Y 


1  =  0-698  b 

2  a  V 


6-28  b 
0)  Y 


where  a  and  00  are  the  real  and  imaginary  parts  of  the  complex  root 
of  stability  equation  (l) .  The  damping  of  the  spiral  mode  of  motion  in 
seconds  is  given  "by  the  equations 


Tl/2 


-0.693  b 

\2  V 


rp  0-693  b 
T2  '  Xg  Y 


ANALYSIS 


Method  for  constructing  a  "boundary  that  defines  a  region  of  satisfactory 
period- damping  relationship.-  For  the  derivation  of  a  method  by  which  the 
period- damping  relationship  of  the  lateral  oscillatory  mode  may  be  expressed 
as  a  function  of  Cnp  and  C the  following  assumptions  are  made: 

1.  The  stability  equation  is  given  in  the  conventional 
form  A>1  +  B\3  +  CA.2  +  d\  +  E  =  0 
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2.  The  coefficients  A,  B,  C,  D,  and  E  are  functions  of  the 
stability  derivatives,  all  of  which  have  prescribed  values  with  the 
exception  of  Cn  and  C7  . 
nP 

3 •  The  stability  equation  has  for  one  of  its  roots 

X-|_  =  a  +  io>  =  Ee*® 


where 


0  =  tan-l  2 
a 


and 


E  =  \[efi  + 


This  root  represents  an  oscillatory  mode  of  motion,  the  period  of 

which  is  given  by  P  =  ^  ^  and  the  time  required  for  the  amplitude  of 
the  oscillation  to  damp  to  one-half  its  original  value  is  given 

ti/2  =  Substituting  Eeie  for  X  in  the  stability  equation, 

and  noting  that  =  cos  n0  +  i  sin  n 0,  results  in  the  expression 

AB^(cos  40  +  i  ain  40)  +  BE3 (cos  30  +  i  sin  30) 

+  CE^cos  20  +  i  sin  20)  +  DB(cos  0  +  i  sin  0)  +  E  =  0  (2) 


Equating  to  zero  real  and  imaginary  parts  in  equation  (2)  gives  the 
equations 

AB^cos  40  +  BE3cos  30  +  CE^cos  20  +  DE  cos  0  +  E  =  0" 

(3) 

AE^sin  40  +  BE3sin  30  +  Cl^sin  20  +  DE  sin  0=0  | 

Multiplying  the  first  equation  of  equations  (3)  by  cos  40,  the  second 
by  sin  40,  and  adding  the  two  resulting  equations  give  the  expression 

A€q  +  +  C^2  +  1^3  +  =  0  (4) 


where 


cos  n0 
E11 


Similarly,  multiplying  the  first  equation  of  equations  (3)  by  sin  40,  the 
second  by  -cos  40,  and  adding  the  two  resulting  equations  give  the 
expression 


A  +  C^2  +  f$3  +  E^4  =  0 


(5) 
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where 


sin  nQ 


Since  0 


are  functions  of  the  terms 


and  a)  of  the  root 


C7  CU.1U.  -L\  CXJ.  X  LUiV/  UXWUW  VJ.  - - - - —  \ 

the  values  of  the  factors  €n  of  equation  (4)  and  of  equation  (5) 

will  he  fixed  for  specified  conditions  of  period  and  damping.  Therefore, 
solving  equations  (4)  and  (5)  simultaneously  results  in  the  0^,0^ 

combination  that  satisfies  any  desired  period  and  damping  relationship 
of  the  oscillatory  mode. 

If  a  suitable  criterion  exists  for  the  damping  of  the  lateral  oscillation 
expressed  as  a  function  of  the  period,  for  example, 

T-L/g  =  f(P) 

this  criterion  will  b©  expressible  in  terms  of  a  and.  ts  by  an  equs.bicn 
of  the  form 

-0.698  b  =  f fax  b\ 


Thus,  for  any  given  period,  there  exists  a  corresponding  value  of  the 
damping  term  a  necessary  to  satisfy  the  criterion.  The  terms  «n 
of  equation  (4)  and  terms  ^  of  equation  (5)  are  expressed  in  terms 

of  R  and  0,  which  are  functions  of  a  and  cn  (r  =  |/a^  +  u> 

and  0  =  tan-1  — \  therefore,  curves  showing  the  variation  of  en 
S'/ 

and  (jf  with  the  period  of  the  oscillation  that  will  satisfy  the  given 
criterion  can  be  constructed.  For  any  value  of  the  period  P  the 
values  of  ^  and  (jf  may  be  taken  from  these  curves  and  substituted 

in  equations  (4)  and  (5)*  Til©  corresponding  C  ,C^  values  may  then 
be  determined  from  a  simultaneous  solution  of  these  equations. 

It  is  important  to  note  that  the  method  described  is  applicable 
to  any  criterion  for  satisfactory  damping  of  the  oscillatory  mode 
expressed  as  a  function  of  the  period.  If  the  criterion  is  expressed 
in  terms  of  the  number  of  cycles  required  to  damp  to  half  amplitude, 
that  is,  V±/2  =  f(P),  the  preceding  method  can  still  be  applied 

Bince  Cl/2  = 

General  procedure  for  use  of  the  method*-  In  the  analysis  of  a 
given  airplane  the  coefficients  A,  B,  C,  D,  and  E  of  equation  (l)  are 
evaluated  in  terms  of  the  known  or  prescribed  values  of  all  the 
stability  derivatives  except  two  ^usually  Cn^  and  which  are  left  in 

symbolic  form.  A  point  on  a  curve,  plotted  as  a  function  of  the  two 
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derivatives  chosen  to  "be  variables ,  defining  the  region  of  satisfactory 
period- damping  relationship  is  obtained  in  the  following  manner:  For  an 
arbitrary  value  of  the  period  P,  values  of  en  and  0n  are  taken  from 
previously  constructed  curves  that  satisfy  the  desired  criterion  for 
period- damping  relationship  of  the  oscillation.  These  values  are  substi¬ 
tuted  in  equations  (4)  and  (5),  together  with  the  expressions  for  the 
coefficients  A,  B,  C,  D,  and  E.  The  two  equations  are  then  solved 
simultaneously  for  the  variable  derivatives .  By  repeating  the  process 
for  successive  values  of  P  the  desired  boundary  will  be  obtained.  In 
order  that  the  significance  of  the  boundary  thus  obtained  be  clearly 
understood  it  is  necessary  to  know  on  which  side  of  this  boundary  there 
exists  a  satisfactory  period- damping  relationship  according  to  the  pre¬ 
scribed  criterion.  A  method  for  obtaining  this  information  is  discussed 
in  appendix  A . 

Application  of  the  method  presented  to  longitudinal- stability  analyses *- 
It  should  be  pointed  out  that  the  methods  presented  herein  are  not  limited 
to  an  analysis  of  lateral  stability  but  can  be  readily  employed  in  the 
longitudinal- stability  analysis  as  well.  If  a  suitable  criterion  exists 
for  the  period- damping  relationship  of  the  longitudinal  oscillatory  mode 
of  motion,  curves  of  and  €n  plotted  against  P  can  be  constructed 

and  points  on  a  boundary  defining  the  region  in  which  this  satisfactory* 
relationship  exists  between  the  period  and  damping  of  the  longitudinal 
oscillation  can  be  calculated.  The  coefficients  A,  B,  C,  D.  and  E  to 
be  used  in  the  simultaneous  solution  of  equations  (4)  and  (5)  for  the 
longitudinal  case  must,  of  course,  be  the  coefficients  of  the  longitudinal- 
stability  biquadratic.  The  variables  of  interest  become,  instead  of  Cnp 
and  any  two  appropriate  parameters  for  longitudinal  stability. 

Method  for  constructing  curves  of  constant  rate  of  spiral  divergence. - 
For  the  combinations  of  Cnp  and  C^p  necessary  to  satisfy  various  cri- 

terions  for  the  period  and  damping  of  the  oscillatory  mode,  undesirable 
spiral  instability  may  be  present.  In  view  of  this  fact  it  is  desirable 
to  know  something  about  the  divergence  gradient  throughout  the  Cnp,C2p 

plane . 

The  usual  method  of  obtaining  this  information  involves  the  solution 
of  a  series  of  stability  biquadratics,  often  a  rather  laborious  procedure. 

The  following  method  is  presented  which  gives  the  desired  information  with 
a  minimum,  of  calculations . 

Assume  that  the  lateral- stability  biquadratic  (equation  (l))  has  as 
one  of  its  roots  X  =  Xg*  This  substitution  for  X  in  equation  (l) 
results  in  the  expression 

Avl  +  E\,23  +  CX-22  +  m2  +  E  =  0  (6) 

where  the  coefficients  A,  B,  C,  D,  and  E  are  functions  of  Cn  and  C}p 
as  before.  If  X 2  =  0,  equation  (6)  is  reduced  to  the  equation  E  =  0, 
which  defines  the  boundary  of  neutral  spiral  stability.  However,  if  X2  is 
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assigned  a  finite  value  representing  a  given  rate  of  divergence,  equa¬ 
tion  (6)  "becomes  a  definite  function  of  Cnp  and  The  solution  of  this 


resulting  equation  for  a  sequence  of  values  of  Cn^  gives  the  corresponding 

values  of  C,  on  a  curve  representing  a  constant  rate  of  divergence  of 

the  spiral  mode  of  motion.  These  curves  can  "be  constructed  for  as  many 
values  of  as  W  he  necessary  for  an  adequate  analysis  of  the 

spiral  instability  of  a  particular  airplane. 


ILLUSTRATIVE  EXAMPLES  AND  DISCUSSION 


In  order  to  meet  the  present  Navy -Air  Force  flying-dualities  requirement 
(references  2  and  3)  the  period  and  damping  of  the  lateral  oscillation 
must  satisfy  the  following  criterion: 


1  /  =  1.5 

1/2 

(0  <  P  <  2) 

> 

1/2  =  2-5F  ■  3-5 

(2  <  P  <  0°) 

—  —  J 

(7) 


This  criterion  is  illustrated  in  figure  1. 


As  an  example  of  the  methods  presented,  calculations  were  made  for  the 
hypothetical  high-speed  airplane  described  in  table  I  by  using  equation  (7) 
as  the  criterion  for  satisfactory  damping  of  the  lateral  oscillation. 

g  » n 

The  values  of  the  factors  - and  - 7TT  ^or  a  sequence  of  values 

(T/t)n  (T/l)1^1 

of  the  period  P  are  presented  in  table  II  and  the  variations  of  these 
factors  with  P  are  shown  in  figure  2.  Note  that  the  ordinates  in  figure  2 

0 

na 

n— 1  n  ' n 


are 


n 


(v/b)n 


and 


(v/b) 


rather  than  e  and  0  in  order  that  the 
n  rn 


curves  in  figure  2  will  apply  for  any  airplane  at  any  flight  condition  if 

the  criterion  described  in  equation  (7)  is  used.  For  a  number  of  values 

of  P.  the  values  of  e  and  0  were  determined  from  these  curves  and 

n  'll 

were  then  substituted  in  equations  (4)  and  (5 ) ,  and  the  corresponding 
C  ,C->  values  were  obtained  from  a  simultaneous  solution  for  these 
nP  zp 
equations. 


The  boundary,  plotted  as  a  function  of  Cn^  and  C^,  defining  the 

region  that  satisfies  the  period- damping  relationship  given  by  equation  (7) 
is  presented  in  figure  3.  For  purpose  of  comparison  the  neutral-oscillatory- 
stability  boundary  (Rj_  =  0)  for  this  airplane  is  also  plotted  in  figure  3* 
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From  inspection  of  the  R^  =  0  "boundary  this  airplane  is  seen  to  he 

oscillatorily  stable  for  almost  any  combination  of  CL  and  C7  shown 

np  ^p 

in  the  figure  0  This  boundary,,  however,  gives  no  indication  of  the  degree 
of  stability  present  at  any  point  in  the  quadrant *  The  curve  defining 
the  Wavy-Air  Force  criterion  for  satisfactory  damping  of  the  lateral 
oscillation  clearly  shows  that  the  criterion  is  satisfied  in  only  a 
relatively  small  part  of  the  stable  region  of  the  CL  ,C7  plane* 

Figure  3,  therefore,  clearly  illustrates  that  the  stability  of  a  given 
airplane  cannot  be  accurately  evaluated  by  merely  constructing  the  R^_  =  0 

boundary  and  noting  on  which  side  of  this  boundary  stability  exists * 

Additional  calculations  were  made  in  order  to  obtain  curves  of 
constant  rate  of  spiral  divergence  for  this  airplane,  and  the  results  are 
presented  in  figure  4.  Should  a  criterion  for  the  rate  of  spiral  divert 
gence  be  established,  it  will  then -be  sufficient  in  a  given  spiraJ.7 
stability  analysis  to  plot  only  the  curve  representing  that  rate  of 
divergence  * 


The  roots  of  the  lateral- stability  biquadratic  usually  indicate  the 
presence  of  one  periodic  and  two  aperiodic  modes  of  motion*  If  this  be 
true,  only  one  cdrve  defining  the  region  of  satisfactory  period- damping 
relationship  of  the  periodic  mode  would  exist*  If,  however,  the  roots 
indicate  the  presence  of  two  periodic  modes  of  motion,  two  or  more  branches 
of  this  curve  may  exist* 


Thus,  reference  4  shows  that  for  an  airplane  equipped  with  an 
automatic  pilot  that  applies  aileron  control  proportional  to  the  dis¬ 


placement  in  roll,  the  roots  of  the  stability  equation  will  represent  two 
oscillatory  modes  of  motion.  Calculations  were  made  for  this  configuration 
by  using  the  stability  derivatives  and  mass  characteristics  presented  in 
table  I  of  reference  4*  Points  on  the  boundary  of  satisfactory  period¬ 
damping  relationship  were  calculated  for  C-,  =  -0*10  by  using  as 


variables  the  directional- stability  derivative  Cn^  and  the  deriva¬ 
tive  C due  to  the  automatic  pilot.  The  equations  and  coefficients 

presented  in  the  section  entitled  "Equations  of  Motion"  were  modified  to 
take  into  account  the  derivative  C.  in  the  same  manner  as  was  done  in 

reference  4* 


Th9  results  of  these  calculations  are  presented  in  figure  5*  For  a 
number  of  points  in  the  Cn  ,C7^  plane,  roots  of  the  stability  equation 


ip 
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were  calculated  in  an  effort  to  understand  more  clearly  the  significance 
of  each  of  the  tranches  of  the  boundary  presented  in  this  figure.  The 
results  of  these  calculations  are  given  in  table  IIIo  At  point  A  the 
roots  represent  two  periodic  modes.,  neither  satisfying  the  Navy-Air  Force 
criterion.  Upon  passing  through  the  branch  of  the  boundary'  in  the  upper 
quadrant  to  point  B*  the  period- damping  relationship  of  one  of  the 
periodic  modes  becomes  satisfactory  while  the  other  remains  unsatisfactory. 
At  point  C  the  relationship  of  the  unsatisfactory  modes  has  improved  some¬ 
what  and  the  other  mode  remains  satisfactory.  At  point  D  the  previously 
satisfactory  periodic  mode  has  become  two  aperiodic  modes*  one  of  which 
is  unstable.  The  period- damping  relationship  of  the,  periodic  mode  is  now 
satisfactory.  At  point  E  the  periodic  mode  has  again  become  unsatisf actory 
and  the  two  aperiodic  modes  are  stable.  Thus*  the  only  part  of 
the  CL  *C7  plane  that  completely  satisfies  the  Navy- Air  Force  criterion 

Up  ip 

is  that  area  on  the  unhatched  side  of  the  branch  of  the  boundary  appearing 
in  the  lower  quadrant.  It  should  be  pointed  out  that  although  the 
oscillatory  mode  satisfies  the  Navy- Air  Force  criterion  within  this  area* 
the  instability  of  one  of  the  aperiodic  modes  is  such  that  the  region 
may  be  of  little  practical  value. 


CONCLUSIONS 


The  following  conclusions  were  made  from  a  theoretical  investigation 
to  develop  a  method  of  calculating  a  stability  boundary  that  defines  a 
region  of  satisfactory  period-damping  relationship  of  the  oscillatory 
mode  of  motion: 

1.  Through  use  of  the  methods  presented  a  boundary  can  be  obtained 
that  defines  a  region  In  which  there  exists  a  satisfactory  relationship 
between  the  period  and  damping  of  the  lateral  oscillatory  mode  of  motion* 
according  to  any  given  criterion  for  this  relationship. 

2.  A  method  is  also  presented  by  which  curves  representing  a  constant 
rate  of  spiral  divergence  may  be  constructed. 

3*  These  methods*  which  are  developed  in  detail  herein  for  the 
analysis  of  lateral  stability*  are  adaptable  as  well  to  the  analysis  of 
longitudinal  stability. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base*  Va.*  February  23*  19^9 
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APPENDIX  A 

A  METHOD  FOR  DETERMINING  ON  WHICH  SIDE  OF  TEE  BOUNDARY 
OF  SATISFACTORY  PERIOD- DAMPING  RELATIONSHIP 
A  GIVEN  CRITERION  IS  SATISFIED 


If  the  stability  biquadratic  has  for  one  of  its  roots  =  a  +  to, 
it  can  be  shown  from  reference  1  that  the  following  parametric  equations 
are  satisfied 


f 


1 


pf3  =  0 


(Al) 


f  -  Hf2  +  ^2f'4  =  0 


(A2) 


where 

2 

[i  =  or 

f  =  Aa^  +  Ba3  +  Ca2  +  Da  +  E 
f±  =  kAa?1  +  3Ba2  +  2Ca  +  D 

f2  =  6Aa2  +  3Ba  +  C 
f3  =  4Aa  +  B 
f4  =  A 

The  coefficients  A,  B,  C,  D,  and  E  appearing  in  the  expressions  for  f, 
fl,  f2 ,  f 3>  and  f4  are  those  of  the  characteristic  stability  equation. 
These  coefficients  are  functions  of  the  stability  derivatives,  two  of 
which  are  assumed  to  be  variable . 

If  the  parameter  p  is  eliminated  between  equations  (Al)  and  (A2), 
the  expression  obtained  is 

(f2f3  -  fif^Jfi  -  f32f  =  0  (A3) 

This  equation  represents  a  curve  of  constant  damping  a,  and  the  frequency 


lb 
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of  the  oscillation  at  any  point  on  the  curve  can  he  calculated  by  using 
equation  (Al)  .  If  a  =  0,  equation  (A3)  reduces  to  the  familiar  expres¬ 
sion  for  the  neutral- oscillatory- stability  boundary  (BC  -  AD)D  -BE  =0- 

Now,  for  every  point  on  the  boundary  of  satisfactory  period  and 
damping,  plotted  as  a  function  of  any  two  arbitrary  derivatives  x  and  y, 
a  root  Xp  =  a  +  to  will  exist  that  exactly  satisfies  the  prescribed 

criterion.  Assume  that  at  a  particular  point  on  the  boundary  (xp,yp)  the 
root  is  M.  =  ap  +  top.  If  a  =  ap  and  x  =  xp  =  xp  +  Ax  be  substi¬ 
tuted  in  equation  (A3),  the  value  7  =  72  >  which  is  located  on  a  curve 
of  constant  damping,  can  be  calculated.  By  substituting  ap,  xg,  and  72 
for  a,  x,  and  y,  respectively,  in  equation  (Al)  the  value  of  n,  and 
hence  o>2,  is  determined  at  the  point  (2*2^2)  ’  if  the  root  ap  +  tog  at 
point  (x2,y2)  satisfies  the  prescribed  criterion,  the  region  in  which 
this  point  is  located  is  the  satisfactory  region  with  respect  to  the 
boundary  defining  the  criterion. 

An  imaginary  value  of  cug  at  point  (xg  ,72 )  indicates  that  no  com¬ 
plex  root  with  the  real  part  equal  to  ap  exists  at  this  point,  Bince  <x>2 
must  be  real  if  the  root  ap  +  to2  is  to  represent  an  oscillation. 

That  is,  ap  +  tog  will  represent  two  real  roots  and  before  the 
point  (xg,yg)  can  be  established  as  satisfactory  or  unsatisfactory,  the 
other. two  roots  of  the  stability  equation  must  be  determined.  However, 
if  Ax  is  chosen  small  enough,  009  will  be  real. 

In  general,  the  satisfactory  region  may  also  be  readily  identified 
if  roots  of  the  stability  biquadratic  are  calculated  at  several  points 
of  interest.  A  method  of  evaluating  the  roots  of  a  quartic  equation  is 
presented  in  appendix  B • 
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APPENDIX  B 

METHODS  FOP  EVALUATING-  THE  POOTS  OF  A  QUAPTIC  EQUATION 


Various  methods  exist  for  evaluating  the  roots  of  high- order  poly¬ 
nomials  (references  5  and  6)  .  Many  of  these  methods,  although  highly 
accurate,  become  rather  laborious  in  actual  application* 

A  method  developed  by  Lin  (reference  5)  and  independently  by  Doris 
Cohen  of  the  NACA  (unpublished)  affords  a  means  of  evaluating  roots  of 
high- order  equations  with  a  minimum  of  computation.  This  method  is 
basically  one  of  synthetic  division  and  under  certain  conditions,  which 
are  discussed  subsequently,  the  desired  roots  can  be  obtained  very 
rapidly . 

In  order  to  illustrate  the  method  consider  the  quartic  equation 


x11  +  BX3  +  cx2  +  DA  +  E  =  0 

As  a  first  approximation,  assume  that  the  equation  has  the  quadratic 

factor  X2  +  jX  +  k  where  J  =  5  and  k  =  -•  Division  may  then  be 

C  C 

performed  as  follows: 


X2  +  (B  -  j)X  +  c  -  k  -  J(B  -  3) 
X2  +  jx  +  k  X^  +  BX3  +  CX2  +  DX  +  E 


[C  -  k  -  J(B  -  J)]X2  +  [D  -  k(B  -  X  +  E 

[c  -  k  -  J(B  -  j^X2  +  J  g  -  k  -  J(B  -  jjx  +  k[c  -  k  -  J(B  -  jT| 


Remainder 


If  the  remainder  is  equal  to  zero  or  negligible,  the  two  quadratics 
obtained  are  good  approximations  to  the  factors  of  the  quartic  equation. 
If  the  remainder  is  not  negligible,  the  procedure  is  repeated  by  using 
as  the  second  approximation  the  factor  X2  +  j  *X  +  k?  where 


l6 
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and 

.  ,  E _ _ 

k  "  C  -  k  -  J(B  -  j) 

This  operation  is  continued  until  the  remainder  is  negligible  •  The  results 
usually  converge  rapidly  if  the  frequencies  of  the  two  modes  are  of  differ¬ 
ent  orders  of  magnitude  or  if  the  quartic  equation  has  small  real  roots . 

The  usefulness  of  the  method  is  greatly  reduced  in  cases  where  the  results 
converge  slowly • 

A  substitute  method  has  been  derived  by  the  authors  to  be  used  as  a 
means  of  obtaining  the  roots  of  the  quartic  equation  if  the  method  of 
reference  5  converges  slowly • 


It  can  be  seen  from  the  division  performed  previously  that 
if  X2  +  jX  +  k  is  a  factor  of  the  quartic  equation,  two  parametric  equa¬ 
tions  must  be  satisfied,  namely 


D  -  k(B  -  j)  =  j(£  -  k  -  j(B  -  jj] 


(Bl) 


E  =  k[c  -  k  -  j(B  -  jj 

The  method  consists  in  solving  these  two  equations  simultaneously  for  j 
and  k.  A  convenient  procedure  is  as  follows: 


Equations  (Bl)  are  rearranged  in  the  form 


k  _  .13  -  £b  +  JC  -  D  = 
2J  -  B 


(B2) 


k2  -  (j2  "  JB  +  C)k  +  E  =  0  J 

The  first  of  equations  (B2)  is  readily  solved  for  k  for  any  given 
values  of  j  .  A  series  of  corresponding  values  of  j  and  k  are  computed 
in  this  way.  For  each  corresponding  j  and  k  the  value  of  the  left- 
hand  side  of  the  second  equation  of  equations  (B2)  is  evaluated-  This 
value  is  plotted  as  ordinate  against  j  as  abscissa.  The  intersection 
of  the  curve  with  the  j-axfs  evidently  determines  the  solution  for  J • 

The  solution  for  k  is  obtained  upon  substitution  of  this  value  of  j 
back  into  the  upper  equation.  The  values  of  j  and  k  thus  obtained 
may  be  substituted  back  into  the  lower  equation  to  check  the  solution. 


The  methods  presented  are  readily  adaptable  to  equations  of  any  even 
order.  In  the  case  where  an  uneven- order  equation  should  appear,  the 
logical  procedure  is  to  obtain  a  real  root  by  synthetic  division  and  there¬ 
by  reduce  the  equation  to  an  even  order. 


NACA  TN  No.  1859 


17 


/ 

REFERENCES 


1.  Brown,  W.  3.:  A  Simple  Method  of  Constructing  Stability  Diagrams. 

B.  &  M.  No.  1905,  British  A. R.C.,  1942. 

2.  Anon:  Specification  for  Flying  Qualities  of  Piloted  Airplanes. 

NAVAER  SR-119B,  Bur.  Aero.,  June  1,  1948. 

3*  Anon:  Flying  Qualities  of  Piloted  Airplanes.  U.S.  Air  Force  Speci¬ 
fication  No.  I815-B,  June  1,  1948. 

4.  Stemfield,  Leonard:  Effect  of  Automatic  Stabilization  on  the 

Lateral  Oscillatory  Stability  of  a  Hypothetical  Airplane  at 
Supersonic  Speeds.  NACA  TN  No.  l8l8,  1949. 

5.  Lin,  Shih-Nge:  A  Method  of  Successive  Approximations  of  Evaluating 

the  Real  and  Complex  Roots  of  Cubic  and  Higher-Order  Equations. 
Jour.  Math,  and  Phys.,  vol.  XX,  no.  3,  Aug.  1941,  pp.  231-242. 

6.  Sharp,  Henry  S.:  A  Comparison  of  Methods  for  Evaluating  the  Complex 

Roots  of  Quartic  Equations.  Jour.  Math,  and  Phys.,  vol.  XX,  no.  3> 
Aug.  1941,  pp.  243-?58  • 


18 


NAG A  TN  No.  1859 


TABLE  I 

STABILITY  DERIVATIVES  AND  MASS  CHARACTERISTICS 
OP  HYPOTHETICAL  AIRPLANE 


W/S,  IB /ft2 
B,  ft 

p,  slugs /ft^ 

Y,  ft/sec  • 

7,  deg  .  . 


^B . 

kXo,  ft  ...  • 

^Z0>  ^  •' •  *  • 
*1  >  deg  •  •  •  * 
Cj  ,  per  radian 

C j  per  radian 
ciip>  per  radian 
Cnr,  per  radian 


Cfp,  per  radian 
Cy  ,  per  radian 
Cyp,  per  radian 


C 

C 

C 

C 


np(fuselage)  ’  »T  radl“ 
Per  radian . 

np(tail) ’  per  radian  *  * 
2^,  per  radian  . 


. 40 

. 20 

.  .  .  .  0.0004 

. 733 

......  0 

. 0.372 

. . 155 

. 2.02 

.  .  .  .  .  9.64 

......  2 

.  .  .  .  -0.197 

.  .  .  .  0.0929 

.  .  .  -0.00732 

~1,kl  Cn3(tail) 

......  0 

.  0 

"1-33  cnp(tail) 
•  •  •  •  •  -0*25 


. Variable 

cnp  "  Cn|3( fuselage) 
.....  Variable 
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TABLE  III 

NATURE  OF  ROOTS  OF  STABILITY  EQUATION 
IN  Cnp,C^  PLANE 


Point 

see  fig.  5) 

Tl/2 

T2 

P 

IO.65 

_ _ 

4-95 

A 

3.08 

“  —  *“  “ 

1.19 

9-95 

_ 

7.05 

B 

3*16 

.  “  “  "  *• 

1-53 

7.72 

_  _  _  _ 

30.2 

C 

3-50 

- - 

2  .16 

1.02 

I.65 

Aperiodic 

D 

4.20 

- - 

3-15 

•793 

.... 

Aperiodic 

E 

2  .50j 

.... 

2.39 

7-93 
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Figure  1.—  Navy-Air  Force  criterion  for  period— damping  relationship  of  the  oscillatory  mode  of  motion. 


Figure  2C—  Variation  of  the  computation  factors 


11 


and 


fii 


with 


(v/*)n 

5  for  the  Navy- Air  Force  damping-period  criterion, 
and. 


(y/b)11- 1 

The  values 


of 


—  are  presented  in  table  II. 


Figure  5,—  Boundary  defining  region  of  satisfactory  period-damping 
relationship  according  to  Navy-Air  Force  criterion.  Airplane  with 
automatic— pilot  coupling  of  aileron  to  angle  of  hank. 


